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An investigation of the spectroscopy of even-even nu- 
clei with random one- and two-body interactions has 
led to the surprising result that the ground state has 
J p = 0+ in typically 60 - 70 % of the cases. This holds 
both for the nuclear shell model |l],|| and the interacting 
boson model ||. To understand this phenomenon it was 
suggested to study the distribution of the lowest eigen- 
values [Q. To this end, a technique was developed in 
|^| for tridiagonal Hamiltonian matrices, such as for the 
one- and two-body interactions in the vibron model. It 
relies on the property that for these matrices it is possi- 
ble to obtain an estimate for the lowest eigenvalue as the 
minimum of a function || . 

The purpose of this comment is to point out that 
for the vibron model the appropriate function is not a 
quadratic polynomial as used in ||, but a more general 
form which in the vibrational basis to leading order in 
l/N, is given by 



az 2 + bz + c + df + ey/{l- z) 2 (z 2 - j 2 ) 



Here z — n p /N and j = J/N ||. In the absence of off- 
diagonal terms (e = 0), the spectrum is that of a three- 
dimensional oscillator characterized by n p = 0, 1, . . . , N 
and J = n p , n p — 2, . . . , 1 (n p odd) or (n p even). In this 
case, the function f(z,j) reduces to a quadratic polyno- 
mial, and the results of || can be used to obtain the 
percentages for ground states with J = 0, J = 1 and 
J = N to be 64.6, 5.2 and 19.8 %, respectively, in close 
agreement with the exact values of 64.4, 4.1 and 24.2 % 
obtained for N — 20 and 100000 runs. In the case of only 
off-diagonal interactions (a = b = c = d = 0), the ground 
state always has J = 0. 

In the general case, both diagonal and off-diagonal 
terms are present. The off-diagonal term introduces a 
deformation in the system, which gives rise to the occur- 
rence of vibrational bands with J — 0, 1, 2, . . . , J max . As 
a consequence, the J = 1 ground states tend to disappear 
in favor of an enhancement of the J — ground states. 
This effect explains qualitatively the change in the J = 0, 
J = 1 and J = N ground states from 64.4, 4.1 and 24.2 
to 72.1, 0.5 and 23.7 %, respectively. 

In order to carry out a more quantitative analysis, 
we have developed a simple method to analyze arbitrary 
functions f(z,j) by sampling the two-dimensional space 
(z,j). As an example, in Figure 1 we show the per- 
centages of ground states with J = as a function of 



x = o\ I ((To + o~\ ) , where <ii denotes the width of the 
Gaussian distributions for the diagonal (i = 0) and off- 
diagonal (i — 1) matrix elements. The case of equal 
widths (discussed in the previous paragraph) corresponds 
to x = 0.5. Figure 1 shows a good agreement between the 
values obtained from an analysis of f(z,j) (dashed line) 
and the exact results (solid line). The discrepancies ob- 
served for the quadratic polynomials of || (dotted line) 
are due to a different treatment of the off-diagonal matrix 
element. 

In conclusion, the approach proposed in || provides 
an interesting method to analyze random interactions in 
terms of random polynomials which is especially suited 
for tridiagonal Hamiltonian matrices. The advantage of 
this method is that it foregoes the diagonalization of 
thousands of matrices. In this comment we have pre- 
sented a careful analysis of the dominance of J p = + 
ground states in the vibron model, and found good agree- 
ment with the exact calculations. 

For cases of interest to nuclear physics, such as the 
IBM or the shell model, the Hamiltonian matrix has a 
more complicated form. Although it can be brought to 
tridiagonal form, the lack of an analytic form for the 
function / leads to complications comparable to those of 
bringing it to diagonal form right away. For these cases 
the method becomes difficult to apply (other than by 
analogy). Many questions remain open to explain the 
robust features observed in randomly interacting many- 
body systems. 
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FIG. 1. Percentage of J p = + ground states as a function 
of x (defined in text) for the exact results (solid), for f(z,j) 
(dashed) and for the quadratic polynomial (dotted) M. 



2 



